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The effective pair potentials between different kinds of dendrimers in solution can be well approxi- 
mated by appropriate Gaussian functions. We find that in binary dendrimer mixtures the range and 
strength of the effective interactions depend strongly upon the specific dendrimer architecture. We 
consider two different types of dendrimer mixtures, employing the Gaussian effective pair potentials, 
to determine the bulk fluid structure and phase behavior. Using a simple mean field density func- 
tional theory (DFT) we find good agreement between theory and simulation results for the bulk fiuid 
structure. Depending on the mixture, we find bulk fiuid-fiuid phase separation (macro-phase sepa- 
ration) or micro-phase separation, i.e., a transition to a state characterized by undamped periodic 
concentration fluctuations. We also determine the inhomogeneous fluid structure for confinement in 
spherical cavities. Again, we find good agreement between the DFT and simulation results. For the 
dendrimer mixture exhibiting micro-phase separation, we observe rather striking pattern formation 
under confinement. 

PACS numbers: 61.25.Hq, 82.70.Dd, 36.20.Ey, 61. 12. Ex 



I. INTRODUCTION 

Dendrimers are highly branched polymeric macro- 
molecules. They are synthesized from monomers that can 
form bonds with at least three other monomers. Starting 
from a suitable central molecule, a number of generations 
of monomers are sequentially added. With the addition 
of each successive generation, the dendrimer becomes in- 
creasingly branched in its internal structureiSiiii^ A 
dendrimer with n generations of monomers is denoted 
a Gn dendrimer. Dendrimers have numerous technolog- 
ical applications - for example in the targeted delivery 
of drugs-'- and, in particular, anti-cancer agents;2i2ii£ in 
light-harvesting applicationsfii and also as a tool for the 
development of nonviral gene deliveryii^ In a dendrimer 
solution, the individual macromolecules feature differ- 
ent conformations due to fluctuations of the monomers. 
However, statistically one finds that the intra-molecular 
monomer distribution is almost spherically symmetric. 
For the statistical behavior of the ensemble of macro- 
molecules only the average over all conformations and 
orientations of the interacting dendrimers has to be taken 
into account, provided no phases with orientational order 
appear. This suggests treating the interaction between 
dendrimers via radially symmetric effective pair poten- 
tials, even though a dendrimer is a complex, structured 
object. Also, given the number of degrees of freedom 
in each dendrimer, determining the structure and phase 
behavior of suspensions of dendrimers starting from an 
atomistic view point is a tough problem. Obtaining an ef- 
fective interaction potential between pairs of dendrimers 
(e.g. between the centres of mass) is a way to surmount 
this. In such a procedure, one is effectively integrating 
in the partition function over the internal monomeric de- 
grees of freedom of each dendrimer, resulting in an effec- 



tive Hamiltonian for the dendrimer mixture that treats 
each dendrimer as a point particleiiii^iiSiiSiii However, 
in principle such an effective Hamiltonian also includes ef- 
fective three-body and higher body effective interactions 
between the particlesii^ii For dendrimers in solution one 
finds that these three and higher body interactions are 
small compared to the two-body potentialSfil so neglect- 
ing these higher body terms in the effective Hamiltonian 
is justified. Given an effective Hamiltonian for point 
particles one can then determine the fluid structure and 
phase behavior by means of theory or simulations. Thus, 
effective interactions are a powerful tool to bridge the 
length scales and to extract macroscopic properties of 
dendrimer solutions from microscopic details. 

In this paper we use such an approach to study binary 
dendrimer mixtures. We find that all effective pair po- 
tentials between dendrimers take a Gaussian formii^ii^ 
The parameters in these potentials depend strongly on 
the specific dendrimer architectures, to the extent that 
we can 'tune' their range and strength by changing the 
number of generations and flexibility between monomers 
in the dendrimers. This enables us to study the occur- 
rence of a number of phenomena that have been pre- 
dicted for Gaussian mixtures, such as bulk fluid-fluid 
demixingi^iSS and micro-phase separationiSi Micro-phase 
separations^ is characterized by a transition in the fluid 
to a state with undamped periodic concentration fluc- 
tuations which is thought to indicate a thermodynamic 
transition to a rather unusual crystalline states^ When 
Gaussian mixtures are subject to geometrical confine- 
ment, pattern formation can spontaneously ariseiSi We 
find dendrimer mixtures exhibiting all these phenom- 
ena. Thus, we are able to demonstrate that the rich 
phenomenology encountered in Gaussian mixtures corre- 
sponds to real systems. 
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This paper is laid out as follows: In Sec. [H] we describe 
our procedure for determining the effective pair poten- 
tials in dendrimer mixtures and give results for two par- 
ticular systems. In Sec. IIIII we briefly describe theories 
that apply for Gaussian mixtures and then in Sec. II VI we 
present theory and simulation results for the bulk struc- 
ture and phase behavior of particular dendrimer mix- 
tures. In Sec. we present theory and simulation results 
for two particular dendrimer mixtures confined in spher- 
ical cavities. Finally, in Sec. lVII we draw our conclusions. 

II. EFFECTIVE PAIR POTENTIALS 

It has been shown by means of monomer-resolved sim- 
ulations and theory, that the effective pair potentials be- 
tween the centres of mass of dendrimers in a good solvent 
can be modelled by a purely repulsive Gaussian potential 
of the form:^^ ^('') — eexp(— r^/R^), where the param- 
eters e > and R depend strongly upon the parameters 
characterizing the dendrimer's specific architecture. In 
particular, they are both affected by the terminal gen- 
eration number and the bond length between successive 
generations. It is therefore possible to 'tune' the effective 
interactions by modifying these parameters. This result 
has also been confirmed by making comparison with ex- 
perimental results»i^ In the present work we determine 
the effective pair potential between different types of den- 
drimers. We find that in binary mixtures, these can also 
be well approximated by Gaussian functions of different 
ranges and strengths. The pair potential between two 
dendrimers of species i and j thus takes the form: 

F.,(r)=e,,exp(-rVR2^-), (1) 

where the parameters eij > and Ry depend strongly 
on the specific properties of the two types of dendrimers 
in question. 

In order to obtain the effective interaction potentials 
between a pair of dendrimers, we employed a simple 
coarse-grained model, the so-called bead-thread modelji^ 
which treats the monomers as hard spheres of diame- 
ter (T, connected by ideal threads of maximum extension 
(1 -I- d)a. The parameter 6 corresponds to the length 
of the polymer chains between two branching points in 
the dendrimer. Increasing the number of chemical bonds 
between two branching points leads to an increase in flex- 
ibility of the polymer and thus to softer but more long 
ranged effective interactions between dendrimers. Large 
values of S are used to model dendrimers with open struc- 
tures and compact dendrimers are modelled by choos- 
ing a smaller value of (5. As the number of generations 
(monomers) is increased, the radius of the molecule, of 
course, increases and the average monomer density in- 
side the molecule grows as well. For a fixed number 
of monomers in the dendrimer, the value of the effec- 
tive pair potential for complete overlap of the centres of 
mass of a pair of dendrimers decreases upon increasing 
6. Upon increasing the generation number G, keeping 




FIG. 1: Effective pair potentials Vij{r) between dendrimers in 
a binary mixture exhibiting macro-phase separation (system 
A). Species 1 is a G4 dendrimers with S = 2.0a and species 
2 is a G3 dendrimers with S — 3. Oct. The corresponding MC 
data are shown as circles. 

5 fixed, the effective interaction becomes more repulsive 
because of the exponential growth of the monomer num- 
ber with G. By varying S and G, one can easily modify 
the dendrimer architecture and thereby systematically 
'tune' the effective pair potential between dendrimers of 
the same kind^ However, it is much less straightfor- 
ward to a priori determine the cross-interaction between 
two dendrimers with different values for the parameters 

6 and G. At the same time, employing mixtures of differ- 
ent dendrimer types dramatically increases the freedom 
to tailor the macroscopic properties of the system. 

To determine effective dendrimer interactions 
Vij{r),i,j = 1,2 in a binary mixture, we employ 
monomer-resolved Monte-Carlo (MC) simulations of 
the bead-thread model for dendrimers of different 
generations and varying thread length S. We simulate 
pairs of dendrimers of species i and j and determine 
in the simulation the probability Pij (r) of finding their 
centers of mass separated by distance r. Thereafter, 
the effective potential Vij (r) can be obtained using the 
relation: 

(3V,,ir) = -\n[g,,ir)], (2) 

where /? = 1 /ksT is the inverse temperature and gij (r) oc 
Pij (r), are the partial radial distribution functionsii^ Due 
to the strong repulsion for complete overlap between den- 
drimers, we use biased MC, as described in Ref. to 
obtain the effective potentials. The MC-results for Vij (r) 
can be very well described by Gaussian functions of the 
form in Eq. (Q, as can be seen in Figs. ^ and |21 

Having derived effective Gaussian pair potentials for a 
given binary dendrimer mixture it is then possible to ap- 
ply theories developed for point particles to calculate the 
mixture structure and thermodynamics. Such a theory 
is described in the next section. Typically, we find that 
most binary dendrimer mixtures exhibit bulk fluid-fluid 
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FIG. 2: Effective pair potentials Vij{r) between dendrimers in 
a binary mixture exhibiting micro-phase separation (system 
B). Species 1 is a G3 dendrimers with 5 = 2.0a and species 
2 is a G4 dendrimers with S — 4.7a. The corresponding MC 
data are shown as circles. 



phase separation (macro-phase separation). In Fig.^we 
display the effective pair potentials for a particular binary 
mixture that does exhibit macro-phase separation. This 
mixture is composed of G4 dendrimers with 6 — 2.0(7 
together with G3 dendrimers with S = 3. Oct. However, 
in one particular mixture that we investigated, we find 
that the mixture exhibits micro-phase separation^ This 
mixture is composed of G3 dendrimers with 5 = 2. Oct to- 
gether with G4 dendrimers with S = 4.7ct. The effective 
pair potentials for this mixture are displayed in Fig. |21 
It can be seen there that the key ingredient for generat- 
ing micro-phase separation is that the effective potentials 
in one species are short-ranged and strongly repulsive, 
whilst the effective interactions in the other species are 
long-ranged with soft repulsion. 



III. THEORY 

Having derived effective pair potentials between the 
different species of dendrimers in a particular binary mix- 
ture, we now use these potentials as input to liquid state 
theories developed for determining the fluid structure 
and thermodynamics. It has been shown that a sim- 
ple mean field DFT which generates the random phase 
approximation (RPA) for the fluid direct pair correlation 
functions^'^ provides an accurate description of the bulk 
and inhomogeneous structure of fluids composed of Gaus- 
sian particles fiSiiSiSS and more generally, for fluids com- 
posed of particles interacting via bounded or even weakly 
diverging soft potentialsiSS*2SiSI The intrinsic Helmholtz 
free energy functional of the inhomogeneous mixture is 
given by: 



where 



drp,(r)[ln(A3p,(r))-l], (4) 



is the (exact) ideal gas contribution to the free energjiS^ 
and Ai is the thermal de Brogiie wavelength for the den- 
drimers of species i. The quantity .Fox[{pi}] is the excess 
(over ideal) contribution to the free energy due to inter- 
actions between the particles. The RPA approximation 
for this quantity is 

^cx[{pi}] ^ ^Yl J ^^'^ J '^^'^ P*(ri)Pj(i"2)Vjj(|ri -ral). 

(5) 

The fluid direct pair correlation functions generated by 
the RPA functional are simply, •^5. 



(2)/ \ 

Cy- (ri,r2) 



-/3y,,(|ri-r2|). (6) 



For a given set of one body external potentials {V^^*'{r)}, 
the fluid one body density profiles are obtained by min- 
imising the grand potential functionaltSiS^ 
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\r))pdr), (7) 



where pi is the chemical potential for species i. 

Given such a DFT, one can calculate the fluid partial 
radial distribution functions gij{r) using the so called 
'test-particle route', in which one fixes a particle of 
species j and then calculates the fluid one-body density 
profiles around this fixed particle - i.e. one sets the exter- 
nal potentials in Eq. Q equal to the fluid pair potentials. 
The partial radial distribution function are then given by 
= where is the bulk fluid density of 

species i. 

An alternative route to determining the bulk fluid 
structure is via the Ornstein-Zernike (OZ) equations.^— 
In Fourier space, the OZ equations for a two-component 
fluid take the form 



D{q) ' 



(8) 



T[{p,}]^T,^[{pi\]+T,^[{pi\] 



(3) 



where hij{q) is the (3 dimensional) Fourier transform 
(FT) of hij (r) = gij (r) — 1 and the numerator functions 
are 

Nii{q) = cii(q) +/?2[ci2('7) - cii((7)c22(g)], 
N22iq) = C22(<z)+/9?[4(g)-cii(g)c22(<z)], (9) 
A^12('Z) = Ci2{q), 

and the denominator 
D{q) - [1 - p^cn(g)][l - P^C22(g)] - Plp^2cl2iq), (10) 
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where c^- (g) is the FT of the bulk direct pair correlation 
function Cij{r)m^ We define the fluid partial structure 
factors as follows: 

S^J{q) ^ S,, + ,/^k,{q). (11) 

In Sec. IIVI we display sonic typical results for the fluid 
partial radial distribution functions and structure factors. 
However, before discussing results for specific systems, 
we remind the reader of some basic connections between 
fiuid structure and phase behavior. 

Given an expression for the Helmholtz free energy func- 
tional for a fluid, it is then straight forward to calculate 
the bulk fluid phase behavior. The Helmholtz free energy 
F for the bulk fluid is simply obtained by substituting 
Pi{r) = p\ into and the pressure P follows as 

P = —{dF/dV)Ni,N2,T, where V is the volume of the 
system. Through a Legendre transformation, the Gibbs 
free energy G{Ni,N2,P,T) = F{Ni,N2,V,T) + PV is 
obtained, along with the intensive Gibbs free energy 
per particle, G/N = g{x,P,T), where x = N2/N and 
N ^ N1+N2. If, at fixed P and T, ^(x, P, T) is a convex 
function of x {g (x) > 0) for all x-values, the homoge- 
neous mixture is stable. If, on the other hand, there are 
parts with g (x) < 0, then the system exhibits bulk fluid- 
fiuid phase separation and the coexisting state points are 
obtained by equating the pressure and chemical poten- 
tials in the coexisting phases. Within mean fleld (van 
der Waals-like) theories such as the present, a fluid- fluid 
coexistence curve (binodal) is always accompanied by a 
spinodal which is the locus in the phase diagram where 
the inverse isothermal compressibility = 0, - i.e. the 
locus of points at which {d^g/dx^)p,T = 0. The com- 
pressibility XT is simply related to the q = value of 
the fluid static structure factor S{q)^ so the spinodal 
can equally be deflned as the locus in the phase diagram 
at which S{q = 0) — > 00. S{q) can be expressed as a 
linear combination of the partial structure factors^^-'^*' in 
Eq. I|ll(l . i.e., as a linear combination of the functions 
hij (q) . From Eq. ^ we see that the spinodal is given by 
the locus in the phase diagram of D{q = 0) = 0. Using 
the RPA approximation aj (r) = —(3Vij (r) in Eq. 1)10(1 , we 
obtain an analytic expression for D[q) and thus we have a 
simple way of determining if a particular dendrimer mix- 
ture has a spinodal and therefore exhibits bulk fluid-fluid 
phase separationtSS 

In addition to exhibiting macro-phase separation, bi- 
nary mixtures of Gaussian particles are capable of ex- 
hibiting micro-phase separationiSi This occurs when at a 
wavenumber q = qc> Q, there is a divergence S{qc) 00 
(and therefore also a divergence in Sij{q) a,t q = qc). This 
is equivalent to D{qc) — > 0. We define the A-line as the 
locus in the phase diagram at which D{q 7^ 0) = 0. This 
line has a U-shape and separates the (p, a;)-plane into 
two regions. Below the A-line D{q) > for all q valueSf^i 
whereas above it the uniform fluid mixture is unstable. 
Using the RPA approximation in Eq. IIUII it is simple 
to locate the A-line for mixtures exhibiting micro-phase 
separation. 



IV. BULK RESULTS 

In this section, we display results for the bulk structure 
and phase behavior of two particular binary mixtures of 
dendrimers, the first of which exhibits bulk fluid-fluid 
phase separation (macro-phase separation) and the sec- 
ond exhibits micro-phase separation. In both these sys- 
tems the macro-phase separation (in the former) and the 
micro-phase separation (in the latter) occur at fluid den- 
sities sufficiently high as to make simulating these sys- 
tems near to phase separation computationally expen- 
sive, particularly since close to the spinodal/A-line, bulk 
simulations become impractical due to the long-range na- 
ture of the fluid correlation functions, requiring very large 
simulation systems to obtain reliable results. Thus, we 
compare results from the RPA DFT with simulations at 
lower densities, in order to evaluate the reliability of the 
theory. Given that the accuracy of the RPA approxi- 
mation improves with increasing density,^'^-^'*'^^ the good 
agreement between theory and simulations that we find 
at lower densities implies that the RPA is reliable for 
determining the phase behavior of Gaussian mixtures at 
higher densities. In our simulation results, we excluded 
finite size effects by performing simulations at the same 
density but different total numbers of dendrimers (20 000 
and 100 000), that both yield identical results. 

A. Macro-phase separation 

The first system we investigate, hereafter referred to as 
'system A', consists of G4 dendrimers with 6 — 2. Oct (in 
the following termed species 1), and G3 dendrimers with 
5 = 3.0(7 (species 2). In this mixture the ranges of the 
interaction potentials Ry are all very similar, but they 
differ markedly in the prefactors e^ . The parameters of 
the effective pair potentials in Eq. for system A are: 

£11 = 11.6 fcsT, Rii = VMcj 
£22 = 3.45 fcsT, R22 = \/3Tcr 
£12 = GAksT, R12 = VSScr. 

The pair potentials with these parameters are plotted in 
Fig. ^ The radial distribution functions gij (r) for sys- 
tem A are displayed in Fig. O for the total bulk den- 
sity pRii = (pi -I- ^2)1^11 = 1-49 and concentration 
^ = P2I P = We plot MC simulation results, results 
obtained by inserting the RPA closure into the OZ equa- 
tions and those from using the RPA DFT via the test- 
particle route. The latter route guarantees gij(r) > and 
is much more reliable than the RPA OZ route. Nonethe- 
less, the RPA OZ results still agree well with the simula- 
tions results for r > Rn. It is at small r where deviations 
and even unphysical negative values of gij{r) appear in 
the RPA OZ results at lower total densities p. 

In Fig.^we display the partial structure factors Sij{q) 
for system A at the state point corresponding to the re- 
sults displayed in Fig. O We see a peak in the struc- 
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FIG. 3: Partial radial distribution functions gijir) for system 
A with a total density pRfi = 1.49 and concentration x — 0.5. 
The circles are MC simulation results, the solid lines DFT 
test-particle results and the dotted lines results from the RPA 
closure to the OZ equations. 



FIG. 5: Phase diagram for the binary dendrimer mixture ex- 
hibiting bulk fluid- fluid phase separation (system A), plotted 
in the total density, pRn, versus concentration, x = f>ilp 
plane. The straight lines are the tielines linking coexisting 
state points on the binodal. 
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FIG. 4: Partial structure factors Sij{q) for system A at con- 
centration X = 0.5 and total density pRfi = 1.49. The circles 
are MC simulation results, the solid lines DFT test-particle 
results and the dotted lines results from using the RPA closure 
to the OZ equations. 



ture factors Sii{q) and £'22(9) at g = 0, indicating that 
the state point for which these results are calculated is 
not too far from the spinodal. System A exhibits bulk 
fluid-fluid phase separation - i.e. we flnd a solution to 
the equation D{q = 0) = 0, the spinodal for this sys- 
tem. In Fig. [S] we display both the spinodal and the 
binodal for system A. These are located at rather high 
(monomer) densities. For example, for x — 0.5 the spin- 
odal is located at a monomer density of a^pm — 0.47, 
where a is the monomer diameter. At such densities the 
reliability of the effective pair potential for dendrimers 
is perhaps questionable. However, the influence of the 
macro-phase separation is observable at lower densities, 



where the structure factors display a local maximum at 
q — Q- see Fig. ^ 

B. Micro-phase separation 

The second Gaussian mixture for which we display re- 
sults, 'system B', is a binary mixture of G3 dendrimers 
with thread length 5 = 2.0(7 (species 1) and G4 den- 
drimers with 5 = 4.7fT (species 2). Here, the 1-1 in- 
teraction is very soft and long ranged, whereas the 2-2 
potential is comparatively strong and short ranged. The 
parameters in the effective pair potentials [Eq. (QJ], ob- 
tained from MC simulations are: 

£11 — Q.iksT, Rii = \/2T(T 
£22 = S.lfcsT, R22 = VlOOcr 
£12 = 3.2 fcsT, R12 = V52(7. 

The pair potentials with these parameters are plotted 
in Fig. 12 The radial distribution functions gij{r) for 
system B are displayed in Fig. for the total bulk den- 
sity pRii = 0.96 and concentration x = 0.5. We plot 
MG simulation results, results obtained by inserting the 
RPA closure into the OZ equations and those from using 
the RPA DFT via the test-particle route. We find good 
agreement between the MG and DFT results. The func- 
tions gij{r) decay rather slowly in an oscillatory manner, 
with wavelength A ~ 2R22. The presence of slowly de- 
caying oscillations indicates that this state point is not 
too far in the phase diagram from the A-line. In fact, 
it turns out that the density /cRfi = 0.96, corresponds, 
a± X = 0.5, to 0.76px{x — 0.5), where p\{x) is the a:- 
dependent value of the total density on the A-line. More- 
over, the fact that 522(0) > 1 indicates that the den- 
drimers of species 2 prefer complete overlap, a further 
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FIG. 6: Partial radial distribution functions gij{r) for system 
B with a total density pRfi = 0.96 and concentration x — 0.5. 
The circles are MC simulation results, the solid lines DFT 
test-particle results and the dotted lines results from the RPA 
closure to the OZ equations. 



FIG. 7: Partial structure factors Sij{q) for system B at con- 
centration X = 0.5 and total density pRfi = 0.96. The circles 
are MC simulation results, the solid lines DFT test-particle 
results and the dotted lines results from using the RPA closure 
to the OZ equations. 



indication of the propensity towards microphase order- 
ing; recall the potential V22{r) is purely repulsive. 

In Fig. 13 we display the partial structure factors Sij{q) 
for system B at the state point corresponding to the re- 
sults displayed in Fig. There is good agreement be- 
tween the MC simulation results and the DFT results. 
We see a peak in the structure factors S'ii(g) and 522(9) 
at g = (7c — 3.5R22^, indicating the propensity towards 
micro-phase ordering. System B exhibits a A-line, i.e., 
we find a solution to the equation D{q = Qc > 0) ^ 0. 
In Fig. IH^a) we display the A-line for system B and in 
Fig-IHtb) the wavelength Ac of the instability, Ac — 2Tr/qc, 
for the points along the A-line. The A-line is located at 
rather high densities: its minimum density pRf^ = 1.23 
corresponds to a monomer density of a'^ Pm — 0.63. The 
length scale of the instability, Ac, is somewhat smaller 
than 2R22 and is rather insensitive to the precise loca- 
tion on the instability line. This indicates that the length 
scale of the A-line instability is mainly determined by 
the larger particles, whose size and interaction range is 
indeed R22. However, it must be emphasized that the 
emergence of this inherent instability is a characteristic 
of the whole mixture and it requires the presence of both 
species with the associated interactions between like and 
unlike species. 



V. CONFINED SYSTEMS 

We determine the inhomogeneous fluid density pro- 
files of systems A and B, from the previous section, 
when these binary dendrimer mixtures are confined in 
a spherical cavity. Much is known about the behavior 
of Gaussian mixtures at planar wall^^ and in spherical 
cavitieS ] i "?'ni? However, this previous work relied (essen- 



6 



4 

Q. 

2 



' 1 ' 1 ' 1 


1 






\ unstable 








stable 




(a) 




0.2 0.4 0.6 
X 


0.8 







T ' 1 ' r 




(b) 

1 1 1 I 1 I 1 I 1 I 1 

■0 0.2 0.4 0.6 0.8 

X 



FIG. 8: (a) Stability 'phase diagram' for the binary dendrimer 
mixture exhibiting micro-phase separation, plotted in the to- 
tal density, pRn, versus concentration, x — p^j p plane. The 
solid line is the A-line. (b) The wavelength Ac = {2n)/qc of 
the instability on the A-line, where qc is the wavenumber for 
which D{qc) = 0. 
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tially) on a guess for the form of the external potentials. 
Here, we determine effective external potentials between 
the dendrimers and the wall explicitly. The potential be- 
tween the individual dendrimer monomers and the cavity 
wall is assumed to be a hard interaction (i.e. the poten- 
tial of the monomer is zero inside the cavity and infi- 
nite outside the cavity). We calculate an effective poten- 
tial between the centre of mass of the dendrimers and 
the wall in a manner similar to that used to derive the 
dendrimer-dendrimer effective pair potentials - i.e. we 
perform biased MC simulations of a single monomer- 
resolved dendrimer near a planar wall where the potential 
of the monomer is zero for z > (z is the Cartesian axis 
perpendicular to the wall) and infinite for z < 0. This 
planar effective wall potential has the following Yukawa 
form: 



iRii cxp[— 2;/Ri] 



for z > 
otherwise. 



(12) 



For mixture A, the parameters are ei = XbQksT, Ri — 
0.34Rii, €2 = rnksT and R2 = 0.39Rii. For mix- 
ture B, the parameters are ei = GSksT, Ri — 0.39Rii, 
£2 = IbOksT and R2 = 0.89Rii. From Eq. ^ we may 
then construct the external potential corresponding to a 
spherical cavity of radius Rwaib 



e.R,iicxp[-(Rwall-r)/R,i] 
(Rwall-r) 

00 



for r < Rwaii 
otherwise. 

(13) 

Of course, this will only provide a reliable approximation 
to the cavity potential when Rwaii ^ ■ 

In Fig. 121 we display the fluid density profiles for sys- 
tem A with Ni = 15000 particles of species 1 and 
N2 = 15000 particles of species 2 confined in a cavity 
with R„aii/R-ii = 12.2. The average total density in the 
cavity corresponds to a state point well inside the fluid- 
fluid demixing binodal. This phase separation is manifest 
in the fluid one body density profiles. There is preferen- 
tial adsorption of species 1 at the walls of the cavity and 
there are almost exclusively just particles of species 2 
for r < 8R11, in the centre of the cavity. Which of the 
two species is preferentially adsorbed by the cavity walls 
is a rather fine balance, determined by an interplay of 
the various external wall potential parameters. Slightly 
changing one of these parameters can lead to adsorption 
of the other species at the wall. 

In Fig. ^1 we display the fluid density profiles for 
system B with A^i = 15000 particles of species 1 and 
N2 = 15000 particles of species 2 confined in a cavity 
with Rwaii/R22 = 7.2. The average total density in the 
cavity corresponds to a state point inside the A-line. As a 
consequence, we see rather striking 'onion-likc' ordering 
in the cavity - alternating layers rich in the two different 
species of particles. The DFT and MC simulation results 
agree well; the DFT is able to capture the details of the 
highly structured profiles, although a slight phase shift in 
the oscillations in the profiles is noticeable. This order- 
ing has its origins in the existence of the inherent A-line 




FIG. 9: Density profiles for system A with TVi = 15 000 parti- 
cles of species 1 and N2 — 15 000 particles of species 2 confined 
in a cavity of radius Rwaii/Rii ~ 12.2. The solid line is the 
DFT result for species 1 and the dashed line species 2. The 
symbols are MC simulation results. 
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FIG. 10; Density profiles for system B with A^i = 15 000 
particles of species 1 and A^2 = 15 000 particles of species 2 
confined in a cavity of radius Rwaii/R22 = 7.2. The solid line 
is the DFT result for species 1 and the dashed line species 2. 
The symbols are MC simulation results. 



instability of the bulk system. Indeed, as can be seen in 
Fig. El the length scale ( that characterizes the onion 
structure is roughly ^ ^ 2R22, precisely as the critical 
wavelength Ac shown in Fig.jSfb). Under planar confine- 
ment, mixture B would have formed alternating layers of 
species 1-rich and species 1-poor lamellae insteadiSi 



VI. CONCLUSIONS 

We have investigated mixtures of dendrimers, finding 
either bulk fluid-fluid demixing or micro-phase separa- 
tion, depending on the choice of the dendrimers' genera- 
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tion numbers and architecture. Here, in contrast to previ- 
ous studies concerning micro-phase separation, we do not 
employ arbitrary Gaussian potentials but those obtained 
from monomer-resolved simulations of dendrimers. Fur- 
thermore, we observe pattern formation under confine- 
ment, finding a strong sensitivity with respect to the form 
of the wall potential. This work demonstrates, therefore, 
that the rich phenomenology encountered in Gaussian 
mixtures can be realized by employing dendrimers as soft 
colloids with tunable repulsions. It would be desirable 
to synthesize suitable mixtures of dendrimers, along the 
lines put forward in thus work, with an inherent A-type 
instability. Under planar confinement, these would orga- 
nize in lamellae, which opens the possibility of using den- 
drimers as lubrication agents, for instance. Another very 
interesting question is the nature of the phase or phases 



that form inside the unstable region enclosed by the A- 
line. A rich variety of one-, two-, or threc-dimensionally 
modulated structures (cyhnders, lamellae, crystals), sim- 
ilar to that encountered in block copolymer blends or in 
ternary oil/ water /surfactant mixtures could be expected 
here as well. We plan to return to this problem in the 
future. 
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